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I.  INTRODUCTION 
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In  a  recent  publication-'-  we  examined  the  susceptibility  of  railgun  arcs 
to  the  Rayleigh-Taylor  interchange  instability.  It  was  concluded  that  these 
plasmas  were  unstable  and  that  the  growth  rates'  were  usually  large  enough  for 
effects  of  the  instability  to  be  seen  in  experiments  of  interest.  In  this 
report  we  will  be  interested  in  extending  the  previous  analysis  to  another 
type  of  Rayleigh-Taylor  effect,  i.e.,  the  current-fi lamentation  instability. 
The  two  types  of  instability  are  differentiated  according  to  the  manner  in 
which  the  initial  steady  state  of  the  plasma  is  disturbed:  In  the  interchange 
case,  the  perturbation  wave  vector  lies  in  the  same  direction  as  the  steady- 
state  current;  in  the  current-f ilamentation  case,  it  lies  along  the  same  di¬ 
rection  as  the  steady-state  induction  field. 


The  basic  model  and  procedure  will  be  the  same  as  for  our  previous  analy¬ 
sis.  In  particular,  we  will  modify  as  needed  the  formalism  worked  out  by 
Tsai,  Liskow, and  Wilcox  ( TLW ) ^  to  treat  both  types  of  instability  in  an  accel¬ 
erating  plasma.  These  authors  considered  a  plasma  in  which  the  conductivity 
varied  linearly  with  the  density,  which  had  infinitely  extended  tails,  and 
which  was  accelerated  electromagnetically  without  a  projectile.  They  concluded 
that  the  plasma  was  subject  to  both  types  of  instability,  with  the  interchange 
case  having  the  larger  growth  rate.  As  pointed  out  in  Ref.  1,  however,  it  is 
more  appropriate  in  railgun  applications  to  assume  constant  conductivity  for 
the  arc.  Under  this  assumption  the  arc  becomes  completely  confined  by  the 
magnetic  field  and  consideration  must  be  given  not  only  to  differences  in  the 
solution  of  the  plasma  equations,  but  also  to  the  boundary  conditions  that 
must  hold  at  the  interface  between  the  plasma  and  the  vacuum.  In  addition, 
both  the  steady-state  and  perturbed  solutions  are  significantly  affected  by 
the  projectile  and  its  presence  must  be  carefully  accounted  for. 

The  basic  model  to  be  used  is  shown  in  Fig.  1  and  is  identical  to  that 
assumed  in  Ref.  1.  The  upper  and  lower  sides  in  the  figure  represent  rails 
which  carry  current  in  the  direction  indicated.  The  current  also  flows  through 
a  plasma  arc  located  in  between  x  =  0  and  x  =  -£.  The  entire  current  configu¬ 
ration  produces  a  magnetic  field  which  interacts  with  the  arc  current  to  pro¬ 
duce,  via  the  standard  force,  a  very  high-pressure  plasma.  This  plasma 
then  accelerates  the  projectile  (shaded  area)  dowmi  the  gun  tube.  It  has  been 
shown 3  previously  that  the  arc  approaches  a  steady  state  in  the  accelerating 
frame  of  reference.  In  this  calculation,  as  in  our  previous  stability  analy¬ 
sis,  the  steady  state  is  subjected  to  perturbations  in  order  to  examine  its 
subsequent  behavior. 


1  Powell,  John  D. ,  "Interchange  Instability  in  Railgun  Arcs,"  Phys.  Rev.  A.  Z4_, 
8262  (1986).  See  also  Ballistic  Research  Laboratory  Report  Ro.  BRL-TR-2276 , 
Jan  1287. 

cIsai,  I.,  Lishor,  D.,  and  Wilcox,  T.,  "Rayleigh-Taylor  Instabilities  of  an 
Accelerating  Inin  Plasma,"  Phys.  Fluids  24 ,  1676  Cl981). 

{  Powell,  John  I.  and  Batter.,  lad  E.,  "Plasma  Dynamics  of  an  Arc-Driver.,  Elec¬ 
tromagnetic,  Probe  czilc  Accelerator,"  J.  Aval,  Phys.  82,  271?  (1881).  See 
also  "P  I  asm  t  1  .nc~ic?  o'  tie  Arc-Driven  hail  Gun."  Ballistic  Research  Labor  a- 
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Figure  1.  Model  for  stability  calculations. 

The  assumptions  indicated  and  discussed  in  Ref.  1  will  also  be  made  in 
this  calculation:  First,  we  will  assume  that  the  rails  are  perfect  conduct- 
crs  and  are  infinitely  extended  in  the  z  direction.  Second,  ve  will  assume 
that  the  temperature  within  the  arc  is  constant  in  both  space  and  time.  The 
spatial  uniformity  has  been  shown  to  be  a  fair  assumption  from  our  previous 
steady-state  analysis,  and  temporal  uniformity  has  been  suggested  as  reason¬ 
able  because  of  the  extreme  susceptibility  of  the  arc  to  radiation.  In  gen¬ 
eral,  temporally  isothermal  conditions  can  be  expected  to  result  if  radiative 
losses  rapidly  smooth  out  temperature  changes  that  would  occur  both  because 
of  Joule  heating  and  because  of  local  expansions  and  compressions  within  the 
arc.  Finally,  and  most  important,  we  neglect  resistive  effects  in  calculat¬ 
ing  the  growth  rates  of  the  instabilities.  In  other  words,  finite  conducti¬ 
vity  is  accounted  for  in  obtaining  the  steady-state  solution  but,  in  the 
perturbation,  the  conductivity  is  assumed  to  be  infinite.  It  is  not  known 
how  reasonable  this  commonly  made  assumption  is  for  railgun  arcs  and  consid¬ 
eration  of  finite-conductivity  effects  is  much  needed.  Some  efforts  under¬ 
way  to  account  for  finite  conductivity  will  be  discussed  in  Sec.  V. 

The  report  is  organized  as  follows.  In  Sec.  II  the  basic  formalism 
needed  to  carry  out  the  calculation  is  presented.  Much  of  this  formalism  is 
also  contained  in  Ref.  2,  but  is  indicated  here  as  well  in  order  to  set  the 
notation.  In  Sec.  Ill  the  calculation  is  carried  through.  In  Sec.  IV  the 
results  are  discussed,  compared  with  the  results  of  TLW,  and  compared  with 
the  results  for  the  Interchange  instability. 
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II.  FORMALISM 


We  choose  the  origin  of  the  coordinate  system  to  be  located  at  the  back 
of  the  projectile  and  to  move  with  it.  All  quantities  calculated  are  those 
measured  by  an  observer  in  the  accelerating  reference  frame.  Let  p,  P  and  v 
be  the  density,  pressure  and  velocity  in  the  arc.  Let  S,  2  and  j  be  the 
electric-field  intensity,  the  magnetic-induction  field,  and  the  current  den¬ 
sity.  Let  c  be  the  conductivity  and  a  the  acceleration.  We  then  assume  that 
the  arc  is  governed  by  the  standard  magnetohydrodynamic  equations: 


+  v  •  (pv)  =  0 

(1) 

>  ->•  ■+•  *-► 
p  —  +  pv  •  Vv  +  V  P  =  J  X  B  -  p  a 

(2) 

■f 

(3) 

-*•  -V  -V  ->■  -V 

V  x  B  =  yj  ,  7  •  B  =  0 

(M 

J  =  c(E  +  v  x  B) 

(5) 

(  1  +  ^  V  A  =  Vr 

r:  t 

(6) 

Equations  (l)  and  (2)  express  mass  and  momentum  conservation;  Eqs.  (3)  and 

are  Maxwell  equations;  Ec.  (5)  is  Ohm's  law;  and  Eq.  (6)  is  the  isothermal 
equation  of  state.  In  Eq.  (6)  a  represents  the  ratio  of  electrons  to  heavy 
particles  in  the  plasma,  kg  is  Boltzmann's  constant,  T  is  the  temperature, 
and  Km  is  the  ion  (or  neutral)  mass.  Physically,  K  represents  the  square  of 
the  acoustical  propagation  velocity  in  the  plasma  under  isothermal  conditions. 
In  Eq.  ' - '  u  is  the  free-space  permeability. 

As  shown  in  Ref.  1  Eqs.  (2)-(6)  can  be  partially  uncoupled  to  produce 
expressions  for  ?  and  v.  One  finds 
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(V  x  B)  +  V2  E' 
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(7) 


v  +  K 


Consequently, 


l:.e  remain ir. 


(I1,  (7)  and  (6)  car.  be  solved  tc  determine 
r.r t  viria'c  16?  ^  ^  ^ y  i  1  f  r^r.  Zcf  •  11  —  ' ' 


According  to  the  standard  procedure  used  in  stability  analysis,  ve  as¬ 
sume  that  the  solution  to  Eq.  (l),  (7)  and  (8)  can  be  represented  by  the 
steady-state  solution  plus  a  small  time-dependent  correction,  i.e.,  for  each 
function  ? 


F{x,t)  =  FQ(x)  +  F^  (x,t) 


(9) 


Furthermore  for  the  current-filamentation  instability,  the  spatial  and  time 
dependence  of  each  is  assumed  to  be  of  the  form 

F1(x,t)  =  F1(x)ei^Z  =  F1(x)  ei(^  .  (10) 


Here  k  =  ka  is  the  perturbation  wave  number  and  ui  the  angular  frequency. 

For  any  given  real  value  of  k,  u  must  be  determined  from  the  calculation. 

If,  then,  Im(w)  >  0,  the  perturbation  grows  in  time  and  the  plasma  is  un¬ 
stable.  The  assumption  that  the  perturbation  wave  number  has  only  a  z  com¬ 
ponent  (lies  along  the  direction  of  the  unperturbed  ^  field)  should  be  con¬ 
trasted  with  the  interchange  instability  treated  previously,  in  -which  k  was 
assumed  to  have  only  a  y  component  (be  along  the  direction  of  the  unperturbed 
current ) . 


The  equations  that  must  be  obeyed  in  the  steady  state  are  obtained  by 
substituting  Eq.  (9)  into  Eqs^  (l),  (7)  and  (8),  retaining  only  the  lowest- 
order  terms,  and  noting  that  v  =  0  in  the  steady  state.  One  finds 


dxa 


d>:  dx 


an 


P  dE 

J-  0 

dx 


rc  ao 


(11) 


(12) 


Similarly,  the  first-order  equations  are  obtained  by  substituting  Eq. 
(9),  with  each  ?  represented  by  Eq.  (10),  into  Eqs.  (l),  (7)  and  (8)  and 

retaining  only  the  first-order  terms.  The  following  equations  then  result: 


■i":l  +  (c0vlx}'  +  1  k  c 0  vlz  =  0 


D  =  —  V  r 
1  x  1>: 


(13) 

(1M 


=  (7.  V 


(15) 


'  \ 
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ij 


0  lx  "1  .  p  lz  U  0  U  o  lx  10 

l“»0vlz  =  lHV7*i  (17) 

where  the  primes  denote  differentiation  with  respect  to  x.  The  expressions 
for  B^  and  v^  are  not  included  since  they  do  not  affect  the  dispersion  re¬ 
lation.  These  variables  can  therefore  be  set  equal  to  zero  without  any  loss 
of  generality. 

Equations  (13)- (17)  can  be  simplified  so  that  all  plasma  variables  are 
expressed  in  terms  of  vpx  by  the  following  operations.  We  first  substitute 
Eq.  (l4)  into  Eq.  (17)  and  the  result  into  Eq.  (13)  to  produce 


ia)pi  =  Z*-=- Fk  lwS(povix}'  + 


k2B 

- 2  BX  v,  ]  . 

U  0  lx 


Similarly,  we  substitue  Eq.  (ll)  into  Eq.  (17)  and  Eq.  (18)  into  the  result 
to  yield 

ikp0viz  =^pr[K(covix>'  *  T  Bo  ''lx1  •  (19) 

Finally,  a  single,  second-order  differential  equation  for  vpx  can  be  obtained 
by  using  Eqs.  (19),  (15)  and  (ll)  in  Eq.  (l6).  The  result  is 


A,v,  +  A,  v:  +  v,  =  0 
1  lx  ^  J.X  o  lx 


■1  uEc  r 


do, 


b-(2.y, 


A  =  —  ( - 
h3  K  v 


and  y  is  give,:  by 


k2B  2  Gnk2  B 

_ 0_  1  .  0 _ o  B' 

UO  us ‘‘-k^K  uKp  0 


Evidently,  then,  once  Eq.  (20)  has  been  solved  for  V]_y,  results  for  Pl»  Vlz’ 
and  B^z  follow  from  Eqs.  (18),  (19),  (lM  and  (15/,  respectively. 

We  next  observe  that  Eqs.  (l)-(6)  hold  not  only  in  the  plasma  but  in  the 
vacuum  as  well  with  p  =  0  and  o  =  0.  They  therefore  imply  certain  boundary 
conditions  that  must  hold  at  the  interface  between  the  plasma  and  the  vacuum 
if,  as  in  the  case  considered  here,  the  plasma  is  completely  confined  by  the 
magnetic  field.  These  boundary  conditions  have  been  derived  for  a  general 
case  by  Kruskal  and  Schwartzschild.  For  the  case  at  hand  we  have 


n  •  [$]  =  0  (23) 

Uj*  =  n  x  [$]  ( 2h ) 

n  x  [?]  =  (n*v) [2]  (25) 

[Ko  +  —■]  =  0  .  (26) 


Here  n  represents  the  unit  normal  to  the  separating  surface  (pointing  toward 
the  plasma)  and  the  brackets  denote  the  change  in  the  quantity  in  question  as 
the  interface  is  crossed.  The  parameter  3*  represents  the  surface  current. 
These  conditions  are  necessary  to  yield  the  value  of  the  derivative  of  v^  at 
x  =  -H0,  needed  to  solve  Eq.  (20)  numerically  as  well  as  to  join  the  vacuum 
fields  to  those  in  the  plasma. 


III.  CALCULATIONS 

In  the  previous  section  the  formalism  needed  to  determine  the  suscepti¬ 
bility  of  the  arc  to  the  current-f ilamentation  instability  was  presented.  In 
this  section  the  calculation  is  actually  carried  through. 

In  Ref.  1  the  solution  of  the  steady-state  equations,  i.e.,  Eqs.  (ll) 
and  (12)  was  obtained  for  c  =  constant.  We  found 


§0  =  ~uj  x/*0  az  ,  -  tQ  <  x  <  0 

=  uj  Sz  >  "  <  x  <  "  *0 

and 

2(1+z)p  -x(l+x/*J 

cn  =  — 7— r  U  d+x)e  0 


(27) 


(2C) 
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In  these  equations,  j  represents  the  current  per  unit  height  on  the  rail  sur¬ 
face,  lQ  is  the  steady-state  length  of  the  arc,  6  is  the  ratio  of  the  arc-to- 

projectile  mass,  and  p^  is  the  area  density  of  the  arc.  The  parameter  x  is 
given  explicitly  by  K 

Uj  B*q  ^  aQZ0 

X  "  2Kp£(l+B)  “  K  »  (29) 

and  must  satisfy  the  transcendental  equation 

1=1%T-  f1  +  X/2  -  ^  (l-e‘x)]  .  (30) 

The  constraint  imposed  by  Eq.  (30)  determines  the  unperturbed  arc  length 

We  now  proceed  to  determine  the  electromagnetic  fields  within  the  vacuum 
to  the  left  of  the  plasma.  These  fields  are  necessary  to  match  the  boundary- 
conditions  implied  by  Eqs.  (23)-(26).  In  the  vacuum  we  have  for  the  first- 
order  fields 

Vx51=0  ,  V  .  S'  =  0  (31) 

■>V 

->  -vV  ^®1  ->  -vV 

7  x  fl  =  *  IT  ’  v  •  h  =  0  •  (32) 

The  superscript  V  has  been  attached  to  differentiate  the  vacuum  fields  from 
those  in  the  plasma.  If  we  make  the  assumption  that  the  spatial  and  time  de¬ 
pendence  of  these  first-order  quantities  is  as  given  in  Eq.  (9),  Eqs.  ( 31 )  and 
(32)  are  easily  solved  to  yield* 

V  -V  k(x+£n)  i d 
B.v  =  B.v  e  0  e 
lx  lx 


"V  k(x+Jt_)  id 
i  B)  e  0  e 
lx 


::  >  C.  It  can  be  sho^i  that  the  resulting 
subse auenilu  in  Fice  2  ana  c,  devend  onlu  . 


disver— 


y  ~V  k(x+2n)  i<t 
E,  =  E  e  0  e 
lx  lx 

V  ~V  k(x+2_)  id 

E*y  =  -  u/k  e  0  e  (3M 

V  ~  k(x+2A)  i(f> 

E,  =iE.  e  0  e 

lz  lx 

~V  'V 

In  these  equations  B^x  and  E^x  are  constants  which  must  be  determined  from 
the  boundary  conditions. 

In  the  steady  state  the  equation  of  the  surface  which  separates  the 
plasma  and  the  vacuum  is  given  by  x  =  -£q  and  the  unit  normal  to  the  surface 
is  simply  n0  =  a^.  In  the  perturbed  configuration,  the  equation  of  the  sur¬ 
face  is 


x  =  -  2  +  6  (z,t) 

0  x 


where  ^  is  given  to  first  order  by 

6  (z,t)  =  — ^  e*^  (36) 

A  Ll' 

and  where  Vq  =  vpx  (x  =  -20).  Consequently,  in  the  perturbed  state  the  unit 
normal  to  the  interface  is,  to  first  order,  given  by 

n  =  a  +  ^  ei6  a  .  (37) 

Xu;  z 

W’e  finally  observe  that  for  zeroth-  and  first-order  quantities  and  we 
have  through  terms  of  first-order 


?,(-*«  +  6  )  =  ?,(- 


1'  "O' 


f  l-i  +  ( - -)  f  V  • 

ov  o  ~x  V  vax  x  =  -i  x 

We  now  use  the  results  of  Eqs.  (36)- (36),  the  plasma  equations  [Eqs. 
(16),  (19),  (1-)  ar.d  (15) j,  and  the  steady-state  results  [Eqs.  (21)  and  (26), 
to  satisfy  the  boundary  conditions  represented  by  Eqs.  (23)~(26).  One  finds 
after  algebra  that  these  conditions  imply  the  following  results: 


6 


where  v '  means  the  derivative  of  v^x 


at  x  =  -  £, 


Evidently  all  constants 


have  now  been  expressed  in  terms  of  v^  and  the  derivative  v',  necessary  to 
solve  Eq.  (20)  has  been  obtained. 


It  remains  now  only  to  solve  Eq.  (20).  In  addition  to  the  condition  on 
v^  at  x  =  -  £q,  it  is  apparent  that  v^x  must  also  satisfy  the  condition 


vx  (x  =  0)  =  C,  (1*0) 

i.e.,  the  plasma  velocity  (in  the  moving  reference  frame)  must  vanish  at  the 
projectile  surface.  For  any  given  value  of  k,  this  condition  can  be  satisfied 
only  for  certain  discrete  values  of  to2  and,  if  any  or  xnese  is  negative,  the 
corresponding  mode  and  hence  the  plasma  is  unstable. 

Our  procedure  for  solving  Eq.  (20)  has  been  to  assume  values  of  k2  and 
vary  until  the  boundary  condition  set  forth  in  Eq.  (ho)  was  satisfied.  By 
continuing  this  procedure  for  various  values  of  k2,  the  dispersion  curve,  in¬ 
versus  k,  can  be  computed.  Prior  to  carrying  out  the  numerical  integration 
we  have  expressed  the  coefficients  in  Eq.  (21)  in  terms  of  dimensionless 
quantities 


-  0- 

Bo  =  uj  ’  Pc  = 


Vcn  —  —  0 

_  k— x  ,  k  =  k£n,  and  u>  =  —=~ 

p  „  ®  JK 


We  have  also  made  the  transformation  z  =  1  +  x/£q  which  maps  the  interval 


(-£q,0)  into  (0,1).  It  is  then  found  that  the  dispersion  curve  for  the  dimen¬ 
sionless  variables,  u.  versus  k,  depends  only  on  6. 


The  numerical  integration  has  been  carried  out  for  a  wide  range  of 
values  of  £.  Typical  of  the  results  obtained  are  those  shown  in  Fig.  2  in 
which  |,.|/y  is  plotted  versus  jk[/x  for  8  =  0.05,  0.1  and  0.5,  and  in  Fi;  .  3 
for  £  =  l.C,  3.0  ar.d  5.C.  Tnese  curves  constitute  the  central  result  of  our 
calculation  ar.I  are  discussed  in  the  following  section. 


Figure  2.  Dispersion  curves  for  p  =  0.06,0.1,  and  0.5.  Respective  values 
and  0.6?. 


IV.  DISCUSSION 


Our  physical  understanding  of  the  dispersion  curves  in  Figs.  2  and  3  can 
be  described  as  follows.  For  the  very-long-vavelength  modes  (small  k),  the 
perturbation  is  hardly  felt  in  local  regions  of  the  plasma  and  the  growth 
rate,  which  is  proportional  to  u>,  decreases  with  increasing  wavelength.  How¬ 
ever,  the  deformation  of  the  plasma  with  it  parallel  to  ^  requires  a  bending 

of  the  lines  of  magnetic  induction  in  order  that  the  S  field  remain  divergence¬ 
less.  The  bending  of  the  lines  is  evident  since  #  0.  The  bent  field  lines 
now  exert  a  stabilizing  force  on  the  plasma  corresponding  to  an  interaction 
of  the  unperturbed  induction  field  with  the  component  of  the  first-order  cur¬ 
rent  produced  by  (see  the  penultimate  term  on  the  right-hand  side  of  Eq. 

( 16 )  ] .  The  strength  of  this  stabilizing  force  increases  as  k  increases  and 
the  lines  become  more  severely  bent.  Thus,  for  sufficiently  large  k,  this 
restoring  force  can  counteract  the  destabilizing  reaction  force  (p^<2q)  and 
the  plasma  is  stabilized.  This  situation  should  be  compared  to  the  interchange 
instability-1-  in  which  the  induction  lines  become  stretched  or  compressed,  but 
not  bent,  and  w  continues  to  grow  as  k  increases. 

The  instability  produced  by  kj|$Q  is  sometimes  referred  to  as  the 
current-filamentation  instability.  This  descriptive  title  is  applicable  be¬ 
cause  the  bent  field  lines  tend  to  "squeeze  off"  or  pinch  areas  of  the  plasma 
in  which  the  current  density  will  be  high  compared  to  that  in  the  bulk  of  the 
gas.  If  the  plasma  is  unstable  the  current  density  in  these  filaments  will 
grow  in  time.  The  filamentary  behavior  is  not  possible  without  a  bending 
of  the  induction  lines  since  the  formation  of  filaments  clearly  requires  ?x$ 
forces  in  both  Cartesian  directions  normal  to  the  direction  of  current  flow. 

It  is  apparent  from  Figs.  2  and  3  that  for  values  of  8  _>  0.5  the  plasma 
becomes  stabilized  at 

k  /'v  C. 5  (Cl) 

cr 

0 . 5 , 


Consequently,  for  fixed  K  (or  temperature)  the  destabilizing  reaction  force 
increases  as  increases  and  larger  values  of  k  are  needed  to  produce  stabil¬ 
ity.  Alternatively,  for  fixed  c„,  the  arc  length  increases  as  T  increases  and 
the  density  at  local  points  in  the  arc  is  reduced.  Correspondingly,  the  re¬ 
action  force  is  also  reduced  and  stabilization  can  occur  for  smaller  values  of 
k.  Tne  sane  result  indicated  by  Eq.  (*-2)  was  found  by  TLW^  for  a  plasma  with 
conductivity  proportional  to  the  density,  which  was  accelerated  without  a  pro¬ 
jectile,  and  which  had  infinitely  extended  tails. 


For  values  of  6  i  0.5,  the  stability  condition  indicated  in  (^2)  becomes 
less  severe  as  can  be  seen  in  Figs.  2  and  3.  Perhaps  the  reason  for  the  re¬ 
laxation  of  the  criterion  at  these  small  values  of  B  is  that  stabilization 
occurs  at  wavelengths  sufficiently  long  that  size  effects  are  important.  In 
other  words,  these  long-wavelength  modes  "see"  that  there  is  a  fixed  boundary 
at  the  projectile  surface  and  the  boundary  exerts  a  stabilizing  effect. 5  in 
general,  the  boundary  becomes  important  for  k  £q-1*  It  is  a  fairly  simple 

exercise  to  show  for  the  curves  in  Figs.  2  and  3  that  kglQ<  1  for  values  of 
8  £  1,  but  not  for  larger  values  of  8. 


It  is  apparent  that  the  maximum  value  of  w  occurs  in  each  of  the  curves 

U3) 


at  points  near  kg/2.  We  may  write  for  these  maxima 
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where  C f  £ )  appears  tc  approach  a  numerical  value  of  about  0.2  as  8  becomes 
large.  That  the  maximum  growth  rate  should  increase  with  a n,  and  decrease 
with  K,  is  evident  from  previous  discussion.  Condition  (Ul*)  should  be  com¬ 
pared  with  the  maximum  growth  rate  obtained  by  TLW,  namely. 


=  C .  2  z  /ft  .  (1»5; 

it  therefore  arrears  that  the  growth  rates  of  the  instability  for  the  constant- 
conductivity  case  are  smaller  than  for  the  case  worked  out  by  TLW,  particu¬ 
larly  for  small  values  of  6.  Undoubtedly,  the  greater  stability  is  at  least 
partly  attributable  tc  the  finite-size  effect  discussed  in  reference  to  Eq. 
f-2'.  A  graph  of  C'£>  versus  £  is  shown  in  Fig.  h  for  the  values  of  8  indi- 
cated  in  Figs.  2  ar.r  5,  as  well  as  for  a  few  additional  ones. 

There  remains  the  question  of  whether  the  instability  is  of  any  import¬ 
ance  in  railgun  experiments.  Evidently,  rapidly  accelerating,  fairly  cold, 
long  arcs  are  the  best  candidates  for  demonstrating  the  effects  of  the  in¬ 
stability,  a  set  of  conditions  not  easily  achieved.  As  examples  for  discus¬ 
sion  we  have  considered  in  Table  I  three  sets  of  data.  The  first  two  corres¬ 
pond  to  actual  experimental  measurements  on  railguns;  the  last  case  is  hypo¬ 
thetical  as  discussed. 


The  first  case  in  Table  I  corresponds  to  an  experiment6  carried  out  by 
Jamison  and  Burden  on  a  small,  one-half-inch,  square-bore  railgun  at  BRL.  A 
2.5  g  projectile  was  accelerated  using  an  aluminum  arc  (Mg  2l  4*6  x  10-26  *<g) 
to  about  a  kilometer  per  second.  Measurements  of  a q,  £q,  the  current  i,  and 

the  muzzle  voltage  Vq  were  made  at  various  points  along  the  gun  tube;  the 
values  indicated  in  the  table  were  obtained  at  a  point  about  15  cm  from  the 
breech.  The  arc  temperature  was  very  roughly  approximated  by  assuming  that 
one-third  of  the  measured  muzzle  voltage  corresponded  to  a  resistive  drop 
across  the  plasma7  and  calculating  the  conductivity  o  of  the  arc.  The  temper¬ 
ature  was  then  inferred  from  the  theoretical  expression^  for  the  Spitzer  con¬ 
ductivity  which  gives  o  as  a  function  of  T.  The  arc  was  assumed  to  be  singly 
ionized.  Once  these  parameters  are  known  K,  x  and  8  follow  from  Hqs.  (6), 

(29)  and  (30),  and  C(6)  can  be  obtained  by  interpolation  using  the  data  in 
Fig.  4.  Finally,  w  and  k  follow  from  Eqs.  (44)  and  (42),  and  the  e-folding 
time  associated  witK  the  mfximum  growth  rate  as  well  as  the  stabilizing  wave¬ 
length  from  the  expressions 


(1*6) 


(1*7) 


7:.e  m:-. Tr.itu::-.  rf  tr.c  last  two  quantities  indicates  that  the  instability  is 
utterly  irit  >  ir.  this  experiment.  Not  only  is  the  e-folding  time  orders 
of  mar r.ituie  larrer  than  the  acceleration  t-'me,  but  also  the  wavelengths  as¬ 
sociated  vitr.  th*  instability  are  large  compared  to  the  railgun  dimensions  in 
the  z  direct  i  :r  (r.zt  accounted  for  ir.  the  model). 
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:.e  second  case  considered  in  Table  I  corresponds  to  diagnostic  data  re- 
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rectlv  from  tr.e  data,  "rut  tr.e  acceleration  was  calculated  from  the  c 
ar.d  (assumed:  ir.iuct&r.ce  grazier.4  of  the  railgun,  namely,  0.1  vdi/m.  °~- 

perature  ar.d  rem ir.ir.,'  j arar.et.ers  were  then  calculated  in  a  manner  identic.-, 
tc  that  for  Cas*:  I.  Note  that  ever,  though  the  arc  is  nearly  an  order  of  mag¬ 
nitude  lonrer  than  before,  the  temperature  smaller,  and  the  acceleration  nearly 
a  factor  cf  tv;  larger,  tr.e  instability  is  apparently  still  unimportant.  In 

ra^ticular,  both  *  ar.d  t  are  still  too  large  for  the  instability  to  have  any 
s  c 

practical  importance. 
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Finally,  ve  consider  a  hypothetical  case,  similar  to  Case  I,  but  for  a 
longer  arc  and  higher  current.  We  choose  the  arc  material  to  be  copper 
(Mq  2l  1-1  x  10-^5  kg),  assume  a  current  of  800  kA,  a  projectile  mass  of  2.5 
g,  and  take  6=1.  The  temperature  was  approximated  by  scaling  the  value  for 
Case  I  according  to  previously  derived  scaling  relations. 9  The  plasma  was 
assumed  to  be  doubly  ionized  here  since  the  temperature  turns  out  to  be  some¬ 
what  higher  than  in  Case  I.  These  parameters  are  then  sufficient  to  specify 
the  steady  state  of  the  arc,  and  the  remaining  variables  in  the  table  can  be 
calculated  as  before.  We  did  use  the  scaling  relations  to  obtain  the  accel¬ 
eration  from  Case  I,  rather  than  use  Eq.  (29),  since  this  theoretical  result 
is  known  to  overestimate  the  acceleration  substantially. 3,9  Obviously,  now, 
the  acceleration  and  arc  length  are  enormous;  indeed,  it  is  doubtful  that 
small  railguns  could  withstand  these  large  currents  and  very  long  railguns 
would  be  necessary  to  accelerate  such  long  arcs.  At  any  rate,  the  e-folding 
time  te  is  sufficiently  short  that  instability  effects  might  be  seen.  How¬ 
ever,  A  is  still  too  long  to  be  allowed  by  any  practical  railgun. 

In  summary,  it  appears  that  the  current-f ilamentation  instability  is  of 
virtually  no  importance  in  railgun  applications,  at  least  whenever  macroscopic 
projectiles  are  accelerated.  Evidently,  the  wavelengths  at  which  stabiliza¬ 
tion  occurs  are  generally  long  compared  to  reasonable  railgun  dimensions,  and 
the  e-fclding  times  are  generally  long  compared  to  typical  acceleration  times. 
This  conclusion  should  be  contrasted  with  that  for  the  interchange  instability, 
in  which  the  growth  rates  always  increase  with  decreasing  wavelength,  and  ef¬ 
fects  of  the  instability  are  always  expected  to  be  important  in  railgun.  app  li¬ 
cations  . 
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Improvements  in  the  present  model  are  not  likely  to  result  in  any  cnan.gf- 
the  basic  conclusion  concerning  the  importance  of  the  current-f ilamentation 
tability  in  railgun  accelerators.  Probably  the  greatest  limitation  of  the 
el  is  the  assumption  of  infinite  conductivity  in  the  dispersion  calcula- 
r.s.  A  prior: ,  however,  one  would  expect  that  finite  conductivity  would 
uee  the  growth  rates  even  more.  This  conclusion  appears  to  be  consistent 
h  recent  were  by  Slcar.^9  and  by  Huerta  and  Decker Sloan  has  carried 
seme  limiting-case  calculations  and  speculated  that  growth  rates  should 
reass  with  decreasin'  conductivity.  Huerta  and  Decker  have  includes  the 
ects  cf  finite  conductivity  in  the  TLW  model  and  worked  out  the  interchange 
e  explicitly.  Although  their  work  is  apparently  still  underway,  early  re¬ 
ts  appear  to  be  consistent  with  Sloan's  speculation. 
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